{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Visualizing event horizon and ergosphere of Kerr black hole\n", "\n", "### Importing required modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import astropy.units as u\n", "import matplotlib.pyplot as plt\n", "from einsteinpy.utils import kerr_utils, schwarzschild_radius" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Defining the black hole charecteristics" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "M = 4e30\n", "scr = schwarzschild_radius(M * u.kg).value\n", "# for nearly maximally rotating black hole\n", "a1 = 0.499999*scr\n", "# for ordinary black hole\n", "a2 = 0.3*scr" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Calculating the ergosphere and event horizon for spherical coordinates" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "ergo1, ergo2, hori1, hori2 = list(), list(), list(), list()\n", "thetas = np.linspace(0, np.pi, 720)\n", "for t in thetas:\n", " ergo1.append(kerr_utils.radius_ergosphere(M, a1, t, \"Spherical\"))\n", " ergo2.append(kerr_utils.radius_ergosphere(M, a2, t, \"Spherical\"))\n", " hori1.append(kerr_utils.event_horizon(M, a1, t, \"Spherical\"))\n", " hori2.append(kerr_utils.event_horizon(M, a2, t, \"Spherical\"))\n", "ergo1, ergo2, hori1, hori2 = np.array(ergo1), np.array(ergo2), np.array(hori1), np.array(hori2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Calculating the X, Y coordinates for plotting" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "Xe1, Ye1 = ergo1[:,0] * np.sin(ergo1[:,1]), ergo1[:,0] * np.cos(ergo1[:,1])\n", "Xh1, Yh1 = hori1[:,0] * np.sin(hori1[:,1]), hori1[:,0] * np.cos(hori1[:,1])\n", "Xe2, Ye2 = ergo2[:,0] * np.sin(ergo2[:,1]), ergo2[:,0] * np.cos(ergo2[:,1])\n", "Xh2, Yh2 = hori2[:,0] * np.sin(hori2[:,1]), hori2[:,0] * np.cos(hori2[:,1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot for maximally rotating black hole" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[,\n", " ]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD4CAYAAAAD6PrjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO2de3Bkd3XnP6dbrcdIo5FGo9E8NJq3H2ODjRHGDlAhQMxAqAzZYlNQW4shbFxJIEsqqUrskF3vAqkKYSsEKjziCmxMlqxxSFi7WBOvIYaEED80BuzxjMejeWhGmpdGGo2s0YxeffaPczXqkfWc7tv39r3nU/Urdf/ubfXp+/r+fuec3+8nqorjOI7jAGSiNsBxHMeJDy4KjuM4zhVcFBzHcZwruCg4juM4V3BRcBzHca5QFbUBxbJmzRrdsmVL1GY4juNUFHv37j2nqq2z6yteFLZs2UJXV1fUZjiO41QUItIzV727jxzHcZwrlEwURCQrIj8Rke8E77eKyNMi0i0i3xSR6qC+JnjfHWzfUvA/7gvqD4rIO0tlm+M4jrM0StlT+DhwoOD9Z4DPqeoO4DzwkaD+I8D5oP5zwX6IyC7g/cBNwG7gSyKSLaF9juM4ziKURBREpB34JeCvgvcCvA34VrDLg8B7g9d7gvcE298e7L8HeEhVx1T1KNAN3F4K+xzHcZylUaqewp8Dvw/kg/ctwJCqTgbve4GNweuNwAmAYPuFYP8r9XN8xnEcxykDRYuCiLwHOKuqe0tgz1K/8x4R6RKRrv7+/nJ9reM4TuIpRU/hTcAvi8gx4CHMbfR5oElEplNe24G+4HUfsAkg2L4KGCisn+MzV6GqD6hqp6p2tra+Ks3WcRzHuUaKFgVVvU9V21V1CxYo/idV/Q/Ak8D7gt3uBh4JXj8avCfY/k9q83c/Crw/yE7aCuwEninWPsdxHGfphDl47Q+Ah0Tk08BPgK8G9V8F/kZEuoFBTEhQ1RdF5GFgPzAJfFRVp0K0zymGfB4uXbJy+TKMjcH4OExMwOQkTE1ZyednPiMCmQxks1BVBbmclZoaK7W1UFdnr534MDEBo6Mz53muc53Pg6oVESvZ7NXnurrayvR5rquzbU6skEpfZKezs1N9RHOI5PPwyitw4YL9feUVGBmxB0RY1042C/X1VlauhMZGK/X14XyfY1y+bOd5eNjO88WLViYmwvvOmho7rw0NM+d51SoXizIgIntVtXN2vR9552rGx2FgAAYHrQwPX93aLwdTU/a9w8Nw6tRMfS4HTU1WWlqgudkfHtdKPg9DQ3aOz5+3MjZWfjumex6Dg1fXNzTA6tVWWlpgxYry25ZS/I5KO6p2Q549a2V4OGqL5mdiAvr7rRw6ZC6KpiZobYW1a+21SNRWxpdXXrFz3N9v53wqxt7ZkRErx4/b+xUrZs5za6v1Jp1QcFFII/m8PRhOnYLTp8N1D4SJ6kwr9+WXzV/d1gbr1tnDI5Pyqb2mBX/6PF+6FLVF187oKPT0WMlkTBg2bLDznctFbV2icFFIE+fPw4kTcPJk5QrBQoyP2+87ccLcSuvWwcaN9gBJUw/i/Hno67PzHIVLKGzyeThzxkomY8KwaZM1BNJ0nkPCRSHpTEzYQ7Knx7rjaWFyEnp7rdTUQHs7dHSYrzqJjI3ZeT5+3ILDaSGft57QqVN2njs6YPNmy2xyrgkXhaQyMgJHjthDMc6+43IwNgaHD1tpaYGtW60XkYRW5cAAHD1q7qEKzyQsmrExizV1d1vvYft2C1Q7y8JFIWkMDdmNcfp01JbEk4EBK3V1Jg4dHZXnk87nzT105Ei8EwOiQtWu/9OnLUNtxw5rBDhLwkUhKQwNwcGDll3iLM6lS7B/vwWot2yBbdviP2huasrcgIcP25gCZ3HOn4dnn7XxD9df7+KwBFwUKp2REThwwHsG18rkpLkbjh41X/SOHfETh8lJOHbMxGB8PGprKpPhYROHpia48UZYsyZqi2KLi0KlMj5uPYOeHvcll4KpKXPH9PSYW2nHjujdSvm82XPoUDKziKJgaAj+7d8s5rBrV3ITD4rARaHSULUHxUsvJTOtNGqmpqzn0NMDO3eaQEQx3uHkSesBjo6W/7vTwJkz5mrdtg2uu85HxhfgR6KSuHABnn/eWjtOuExMWMyhp8dalOXyRQ8Nwb595gt3wkXVXHInT8LNN3u8IcBFoRLI5y0g2t3trqJyc/Gi+aJbW+E1rwlvUr7xcesZTE/r4JSPS5fsHG/YYOe4ujpqiyLFRSHuXLgAP/mJzVvjREd/P/zgB+ZS2rGjtC6lEyesV+JB5Gg5eRLOnYNbbkl1r8FFIc4cPmyxg3LPUurMTT5vwf2+Prj1VsuBL4bRUfjZz+xB5MSD8XHrNXR0mEsphRPvuSjEkbEx6x34+tPxZGQE/vVfLQh9443L7zVMJwvs3++jzePK8eM2meDrX29jHFKEi0LcGByEvXt9cFLcUbUU1rNn4bbbbGGYpXD5sgm+9w7iz8gI/OhH8NrX2txZKaFox6iI1IrIMyLyMxF5UUT+e1C/VUSeFpFuEfmmiFQH9TXB++5g+5aC/3VfUH9QRN5ZrG0Vx7FjlkPtglA5TD84lpIEcOqUxSVcECqHqSkT8X37UuPGLUW0bAx4m6reAtwK7BaRO4DPAJ9T1R3AeeAjwf4fAc4H9Z8L9kNEdmHrNd8E7Aa+JCLpcOipWqrpCy+k5sJLFPm8ZQ499dTcg8ympix20NXlY0sqlaNH4emnU5EMULQoqDE9J3MuKAq8DfhWUP8g8N7g9Z7gPcH2t4uIBPUPqeqYqh4FuoHbi7Uv9kxM2MXW0xO1JU6xnDsHP/zh1bGgkRH4l3/xVNMkcO6c9QoTPjV5SWIKQYt+L7AD+CJwGBhS1clgl15gY/B6I3ACQFUnReQC0BLUP1Xwbws/M/v77gHuAejo6CjFT4iGy5dNEHymy4omn58pOjaG/vBpprbugBUryB7Yh+SnELGZujOZVCa0JIeLF00Ybr+9+OyzmFISUVDVKeBWEWkCvg3cUIr/u8D3PQA8ANDZ2VmZo7lGRszdUMlLJCaYfN48BWNjM38nJmBiEiYn7PXUlJVXhxIUfnaIKakiq3NkFwlkMzazQjZrUyxNl+pqqK6Bmmqbl88FJIaMj1vsr7PTVntLGCXNPlLVIRF5ErgTaBKRqqC30A70Bbv1AZuAXhGpAlYBAwX10xR+JlkMD9tFlQL/ZNyZmjJdvnQJRi/BpVETgLFxzAkaBjojKIuRy0FtrZW6FbCiztaw96l6ImZqysYz3HYbrF8ftTUlpehLS0RagYlAEOqAX8SCx08C7wMeAu4GHgk+8mjw/t+C7f+kqioijwJ/KyJ/BmwAdgLPFGtf7Bgash6CBxzLjqqNFxsZsTI6CpcuE97DvwRMBL2S2QPaq6ttxo36eli50oTCexVlJp+39PHXvc7WAk8IpWhvrAceDOIKGeBhVf2OiOwHHhKRTwM/Ab4a7P9V4G9EpBsYxDKOUNUXReRhYD8wCXw0cEslBxeEsqJqLuDhYSsXLyZnrNj4uJXpefNETCAaG00kVq6MZnLX1KEKzz1nrxMiDKIVPsFaZ2endnV1RW3G4gwPw49/7IIQMhMTNl3U0JD9jVIELKYwufiOISBiwtDUDM1N8Vs3KHGImCtpw4aoLVkyIrJXVTtn17tnshxcvOg9hBAZG7MW8+AgjFwk1u6gcqE600M63mNLUjc32zr2K1ZEbV0Cme4xVFVVfPDZRSFsLl+2oLKvnFVSJiZMBAYGLD7gLMx0MP3kSQtat7TYipTegyghqjZA8c47Kzpd1UUhTCYnbRyCp52WBFUYumBjwy4M+dIS18rlyzbRa1+frUbZ2mo9CA9Ul4CpKbvn3/zmil3q00UhLKZbDT4wrWgmJmzeubNn3QNXaqYzsY4ft97D2jZLe3WKYHqWgre8pSIX7HFRCIt9+3zq6yK5eBFOnzY3kfcKwmVqakZ4GxuhbR00rbL4qXMNjI7aOIY776y4NDAXhTDo6bEZT51r4sIFm1DUO1nRMB2grquzBchaWiruuRYPBgdtkstbbonakmXholBqBgetl+Asm/NDcLIv8fONVQyXLtnkoH19lmm5Zo2Lw7I5fhyammDz5qgtWTIuCqVkbMxGOPr018viwgXo7XUxiCvj49bxPXnSxmetWeNupWWxb58twtTUFLUlS8J1v1So2mIcvkDOkhkdtSWoDx50QagExset57Bvn/XqnCWSz1fUWhouCqWiu9sDy0tkYiJ4uLzocYNK5NIlOPSyCfqoZ1svjUuXbKGlCsBFoRQMDVlz11kQVTh9xhaZ6+/HRx5XOMPD1mvoOZ6cOaVC5dSpilhsyWMKxTI1ZcPbPWdyQS5etN7B6GjUljglReHMaRgcgC1bKnogb3nYt8/Suerro7ZkXrynUCwHDrhDfAHyeThxAl7c74KQZCYm4NAhKxXiOo+GqSn46U9j3Yh0USiGwUFr/jpzMjJiDaNTp3BXUUo4f95S8wcHo7YkxsT8ueGicK3k8xUTOCo3qpa+uP+AJ2OlkclJy7s4csRjDfPy0kux7Tq7KFwrhw759JxzMD5uMffeXrx3kHLOnbOeontX52BqyrpUMcRF4Vq4eNGaQs5VDA/Di55m6hQwNgb799ucSs4szp61yb1iRtGiICKbRORJEdkvIi+KyMeD+tUi8oSIHAr+Ngf1IiJfEJFuEXleRG4r+F93B/sfEpG7i7UtNPbt81HLszh9Gl466EFG59Wo2ojoo0f9tnkV+/bFzsdWip7CJPB7qroLuAP4qIjsAu4Fvq+qO4HvB+8B3gXsDMo9wJfBRAS4H3gjcDtw/7SQxIozZ7zZU0A+D4cPB+nX7i5yFqC/31zp3nAo4NIlu4FiRNGioKqnVPW54PUrwAFgI7AHeDDY7UHgvcHrPcDX1XgKaBKR9cA7gSdUdVBVzwNPALuLta+kqFpf2AHs5n7pJVv9zHGWwsiIuRh9JHQB3d2xysgoaUxBRLYArwOeBtpU9VSw6TTQFrzeCJwo+FhvUDdf/Vzfc4+IdIlIV385p5Y4ftyDywHTvmI/HM5yGR+HA/s99nSFqalYzYhQMlEQkQbg74HfUdWrTreqKiV0LqjqA6raqaqdra2tpfq3CxOzExclo6MmCL7stHOtTN9OPp4h4MSJ2LSwSiIKIpLDBOEbqvoPQfWZwC1E8HfaEd8HbCr4eHtQN199PDh61J+C2HV74ID7hZ3iUYXuw5a6mnpUzRcbA0qRfSTAV4EDqvpnBZseBaYziO4GHimo/2CQhXQHcCFwMz0O3CUizUGA+a6gLnomJ2MXDIqCV16x1l3MkiWcSkZtkJtPMExslhssxYR4bwL+I/CCiPw0qPtD4E+Ah0XkI0AP8KvBtseAdwPdwCjwYQBVHRSRTwHPBvt9UlXj0bk8dswcoSlmZAReftkFwQmH6VkfyuUNji0vvwydnZGaULQoqOqPgPnWYXr7HPsr8NF5/tfXgK8Va1NJmZqypkyKGR31HoITPkeP2nKfLS1RWxIhp05ZC6yhITITfETzYpw4kepYwuXLLghO+Th8xJZnTTURu6pdFBZCNdW9hIkJEwQPKjtlQ21asVTPl9TbG+m4BReFhThzJrVXZz5vN2eKO0lOROTz5lpP7bWXz1scMyJcFBYixnOeh83Ro7FJm3ZSyPSiPal1Wx4/HtlEUS4K8zEyktoE6pMnfeoKJ3pGR1PcLhsbC1anKj8uCvPR0xO1BZEwPAy98Rky6KScwcFYzi5dHiJ6BrkozEU+H6wSky4mJoLEB5/t1IkRx+MzA0R5GRiI5Ie7KMzFmTOpHKx2+LBnGjkxRO3anJyM2pAIOHFi8X1KjIvCXERwIqLm9JlYjLB3nDkZG0upR7e311Ljy4iLwmzGx1M3Ecvly9CbPh10KoyBATh/Pmoryszly2XP+nBRmM2pU6lbM/DIkdT9ZKdCOXYshW6kvvJmfrgozCaiNLCo6O9PaRDPqUgmJlLo3T19uqwuJBeFQiYmUjU2YWIiWFvZcSqI1DVkxsfL6kJyUSjkzJmyB3WipK8vxSNGnYqmpydVt2pZB2u4KBSSolEyo5fgbLri6U6CuHgxZUt5njlTtq9yUZhGNVVZR70n8EFqTkVz4kSKEiRGR8vmMyvVGs1fE5GzIrKvoG61iDwhIoeCv81BvYjIF0SkW0SeF5HbCj5zd7D/IRG5e67vCo3z51OT1jAyAkNDUVvhOMWRuuzxs2cX36cElKqn8NfA7ll19wLfV9WdwPeD9wDvAnYG5R7gy2AiAtwPvBG4Hbh/WkjKQoqurjJnuDlOaJw8maLeQpmeUSURBVX9Z2C2h28P8GDw+kHgvQX1X1fjKaBJRNYD7wSeUNVBVT0PPMGrhSY8UpJ1dPGir2zlJIdUJQwODJRFAcOMKbSp6nTS/2mgLXi9ESjMNO4N6uarfxUico+IdIlIV38p1HNqKjX+lJQNw3BSwKlTKclEmpoqS4uuLIFmVVVKGNZU1QdUtVNVO1tbW4v/h+fPp6IPOjYGg2mbJsBJPGNjMJSW3m8ZxiuEKQpnArcQwd/pKEkfsKlgv/agbr768EnJhCr9/XjGkZNIzqQlm7wMz6owReFRYDqD6G7gkYL6DwZZSHcAFwI30+PAXSLSHASY7wrqwicFopCyjFsnZQwPp2RN5zI8q6pK8U9E5H8DbwXWiEgvlkX0J8DDIvIRoAf41WD3x4B3A93AKPBhAFUdFJFPAc8G+31SVcszPCUF8YShC75WgpNs+vuhvT1qK0JmbAwuXYK6utC+oiSioKofmGfT2+fYV4GPzvN/vgZ8rRQ2LZnLl1PRxBhIS4aGk1oGBlIgCmDB5hBFwUc0pyA/M0XJVU6KGRtLyUR5IT+zXBRSsNzY0FAqkqscJx3zIYX8zHJReOWVqC0IHe8lOGkhBTkjoT+zXBQSLgqqLgpOehgbszBhohkdDbXr76Jw8WLUFoTKxYu+ZoKTLhIfJlQ1YQiJdIvC5cuJf2KmIGTiOFeRims+xMZsukUhRLWNC6nIxnCcAlJxzXtPISQuXYragtBJxQ3iOAVMTKRg6JGLQkgkXBTGxlKzbpDjXEXCQ4WhRtPTLQoJb06kwDvmOHOS+GvfRSEkEp67lvCOkOPMy2jSr/0QG7TpFoXx8agtCBUXBSetXE76te+iEBIJF4WEe8ccZ14Sf+1PToY2gC3dopDwuaQTrnmOMy+qib+9Q8siSbcoJDw1ZzzpN4XjLEDiRSGkH+iikFAmJ/GlN51Uk3hRCGk2hvSKQj5vfcyEkmC9c5wlkfh7IC3uIxHZLSIHRaRbRO4N7YsSvsBAwn+e4yxK4u+BNASaRSQLfBF4F7AL+ICI7ArlyxLcS4AU3BCOswhTSb8HQnqGxUoUgNuBblU9oqrjwEPAnlC+KeGikPCf5ziLk/R7ICWisBE4UfC+N6i7ChG5R0S6RKSrv7+/bMZVEi4KjuNcC3EThSWhqg+oaqeqdra2tl7bPxEprVExI+E/z3EWJ+n3QEg3edxEoQ/YVPC+PagrPQl/ambidmYdp8xkkn2Lp0YUngV2ishWEakG3g88Gso3Jfypmc1GbYHjREvi74GQnmFVofzXa0RVJ0XkY8DjQBb4mqq+GMqXZTKmtAl1vif+hnCcRUj8PRDSD4yVKACo6mPAY2X5smw2sSNcqqown2oyNc9xFqUqdk+3EhPSD0y2D2UxEnzViEAuuT/PcRalujpqC0LGRSEEcrmoLQiVxN8UjjMfkvjbO7QfmG5RSPhTs7Y2agscJxpqqhOeYJjNhhZodlFIMDU1UVvgONGQ+AZRiM+udItCwp+adXVRW+A40ZD4az9E1Uu3KCS8ObFiRdQWOE40JP7aD7FBm25RSHhzorY2BbnajjMH9fVRWxAyIT67XBQSjEgKbg7HmUU2m3gngItCaCS+jwkNDVFb4DjlpaEh4ZlHEOqzK92iUFub+DmQGhujtsBxyksqrvkQXQDJfiIuRgr8Kw0Nidc9x7mKVauitqAMuCiESML9K5lMSlpOjoMN8k28V7iuLtQMEheFlSujtiB0mpujtsBxykNTU9QWlIGQn1kuCiloRjc1kfxVqBwHWL06agvKQMjPLBeFFIhCLgeNye8QOSknl0vF7eyiEDr19YmeQnuaNWuitsBxwmX16hSkokLoPjIXBUiFI7K52Uc3O8mmtTVqC8pAGSLpRYmCiPx7EXlRRPIi0jlr230i0i0iB0XknQX1u4O6bhG5t6B+q4g8HdR/M1ijuTykIBKbzUJLS9RWOE441NenIOsIrAEbcneo2J7CPuDfAf9cWCkiu4D3AzcBu4EviUhWRLLAF4F3AbuADwT7AnwG+Jyq7gDOAx8p0ralkwJRAFjbFrUFjhMOa9dGbUGZKEMkvShRUNUDqnpwjk17gIdUdUxVjwLdwO1B6VbVI6o6DjwE7BERAd4GfCv4/IPAe4uxbVmkxBm5oi4lgTgnVeRyKeoFx10UFmAjcKLgfW9QN199CzCkqpOz6udERO4RkS4R6erv7y/e2lwuFeMVANavj9oCxyktbW0pGbWfyZTFq7HooRSR74nIvjnKntCtmwdVfUBVO1W1s7VU0aWUpOesWpX4mT2cFJHNpsh1VKZskUVzMVX1Hdfwf/uATQXv24M65qkfAJpEpCroLRTuXx5aW+HIkbJ+ZVRs2AiHXo7aCscpnra2VGSUG2VKrwqr0/Uo8H4RqRGRrcBO4BngWWBnkGlUjQWjH1VVBZ4E3hd8/m7gkZBsm5uWlpT0QaG5KfFTPjkpIJuFdeuitqKMVIIoiMiviEgvcCfwf0XkcQBVfRF4GNgP/CPwUVWdCnoBHwMeBw4ADwf7AvwB8Lsi0o3FGL5ajG3LJptNjQsJoL09agscpzg2bkxRL6GmpmzTvxZ1SFX128C359n2x8Afz1H/GPDYHPVHsOyk6Ghrg7NnIzWhXDQ2QlMzDJ2P2hLHWT41NSmKJYD92DJlSKbDX7JUUtUXhc0dqcjEdRJIx+bUeHuNMj6b0nRYF6e2NjUD2cBaWxs2RG2F4yyPpiaLi6WGbLasc3i4KMwmZYn869enYJFzJzFks7B5c9RWlJm2trJOXOaiMJuUNZ0zGdi2DV9vwakI2tuth5sqNs47jjcUXBRmU1eXojHzRkMDrE9XOMWpQBobUxZcBpttocw/2kVhLlKYr7lxY0pmmXQqkmwWtm5NYWLEhg1lj6i7KMzFhg2pW3wgk4Ht21OW0eFUDNu2pdBtBNDRUfav9EfAXFRVpS62AOY527Ilaisc52ra2lKVFDhDY2MkC4C5KMxH6lIcjDVrUui3dWJLQwNs2rT4fokkomeQi8J8NDeXbVh53Ojo8LmRnOjJ5WDHjpS6NKuqIottpvFwL52tW6O2IBIyGdi5M6U+XCcWZDJw3XVQXb5FeeNFe3tkEzu5KCzExo2pfTLmcnZTpize7sQBsaSHVK/7sW1bZF/torAQmUyqI691dSYMqey+O5GxZXNKA8vTrFsXqSL67b4YW7emaH7eV7NyJWzfkcL8cCcS2ts90YEdOyL9eheFxcjlUpuJNE1zk0+F4YTP+vWpzAS/mjVrIu8muSgshe3bU+9cb2mB7S4MTkisW5fi1NNCdu6M2oKiV177rIi8JCLPi8i3RaSpYNt9ItItIgdF5J0F9buDum4RubegfquIPB3UfzNYrjMe1NSkOrYwzbQwuCvJKSUbNkQycDd+tLTEYvXHYnsKTwA3q+prgZeB+wBEZBe2/vJNwG7gSyKSFZEs8EXgXcAu4APBvgCfAT6nqjuA88BHirSttOzYkerYwjQtLbBjpwefndLQ3p7Kqcbm5oYborYAKFIUVPX/BesuAzwFTJ/ePcBDqjqmqkeBbmypzduBblU9oqrjwEPAHhER4G3At4LPPwi8txjbSk51tbmRHJqb4PrrU+9Rc4pBLIcj9TGEadraYPXqqK0AShtT+DXgu8HrjcCJgm29Qd189S3AUIHATNfPiYjcIyJdItLV399fIvOXwPbtviJNwMqVsGtXaodxOEWQzcJ1O8u6mFi8EYEbb4zaiissKgoi8j0R2TdH2VOwzyeASeAbYRo7jao+oKqdqtrZWs4rK5u1JrID2DiGXbt8Sgxn6VRX2/Mvgnne4svmzdbKigmLOslV9R0LbReRDwHvAd6uqhpU9wGFuQTtQR3z1A8ATSJSFfQWCvePF5s2wbFjcOFC1JbEglzOXKE9PVDOTptTeaxcaaG5XC5qS2JELhe7hmax2Ue7gd8HfllVRws2PQq8X0RqRGQrsBN4BngW2BlkGlVjwehHAzF5Enhf8Pm7gUeKsS00ROA1r4nailiRyZh/eMsWz0xy5qZtnT37XBBmccMNsZvgqdiYwl8AK4EnROSnIvIVAFV9EXgY2A/8I/BRVZ0KegEfAx4HDgAPB/sC/AHwuyLSjcUYvlqkbeHR3Ow5dHOwdq3HGZyryWatd7C5wzPWXsWqVbEcGCszHp/KpLOzU7u6usr/xRMT8OSTMDZW/u+OOVNT5mEbGIjakmiZkiqyV3In0kdDg+VmeCNhDkTgLW+JdHp+Edmrqp2z6127r5VcDm6+OWorYkk2aw8DHwieUsTGHtx4owvCvGzbFtv1Wnw0VjFs2AAnT8KpU1FbEktaWiy4eOwYDA1FbY1TDlassOfdihVRWxJjGhpiF1wuxHsKxfLa13pzaAGqq2367e3bPciYZCToHeza5YKwICJw662x7kK7KBRLdbUJg7MgLS2WtNXaik+qlzBWrbJzu2GDB5MXZefOyGdBXQw/haVg3TqfMG8JVFVZ6upNPuAtEdTU2DPu+ut9oP+SaG6OxSyoi+ExhVJx000wOAjDw1FbEnvq683NMDAAvb2ewFVpZLPWK2hr857Bksnl4LbbKuKAxd/CSiGTgde/3mdSXQbTLqWOzR5vqAQyGVsI55Zb7G8FPN/iw623VkywxZ9gpaShwe6YvXujtqRiyGRgXRu0roGzZ+H0aRsC4sQHEVjbBhvWu3hfEzt2mIu5QnBRKDUbNti8SN3dUVtSUWSz1vpcu9bmUDp9GsbHo7Yq3WSzdj7WrXMxuDMs+AIAAA+dSURBVGZaW2OzTsJScVEIgxtugFdegTNnorak4shm7SG0dq2FaE6fhtHRxT/nlI7qaosXtLa6N7Qo6uvNpVxhE4L5KQ8DEQsq/ehHJg7OsslkbGXCNWvsEJ49ayJR4bOyxJrGRhPj5uaKe47Fj1wO3vjGiuxiuSiERVWVXRQ/+hFcvhy1NRXNypVWOjrg3DlzL/khLQ25nAX81671tNKSkcnAG95gPYUKxEUhTOrqTBh+/GOPnpaAXM7iDuvXw8iIpbQODMBkeuecuyYyGesNtLTYwDPvFZSY173ODm6F4qIQNo2N1mp46inI56O2JjE0NFjp6LChIYODNr+Sa+/cZDImAKtX26pnMZ5lobKZHtpdwbgolIOWFhOGZ591YSgxIvawW7XK4g0jIyYOFy54gLq62o5LUzOsavRxBaFz442JmNnARaFcrF1rmQhdXR4tDQmRmfjDpk2W0jo8PFOSnuKazcLKRuucrmo076VTJq67zsYjJAAXhXKybp0Jw3PPeY+hDFRXz2QwgU2nMTIyU0ZHK1uf6+osltnQAA0roa7W4wORcN11sZ4Ke7kUJQoi8ilgD5AHzgIfUtWTIiLA54F3A6NB/XPBZ+4G/ij4F59W1QeD+tcDfw3UAY8BH9dKXxZuLtavh85O6zG4MJSVmhor0zHAfN6ymC5ehNFLcGnU3setR1FVZQJQW2szJUwXjwvEgBtuqIhJ7pZDsT2Fz6rqfwEQkf8M/FfgN4B3ATuD8kbgy8AbRWQ1cD/QCSiwV0QeVdXzwT6/DjyNicJu4LtF2hdP2tosK+nZZz11JkIymZkHbCFTUzA2DmOXrXcxNg4T4yYWExN2yqamSmdDLmelKgfVOevhTAtYba0PIIstr3lNImIIsynqclPVwilB67EHPVjv4etBS/8pEWkSkfXAW4EnVHUQQESeAHaLyA+ARlV9Kqj/OvBekioKYD6Nn/s5ePppnyY0ZmSzsKLOynyomjjk8yYQU1NWl8/bTTDVuh6tW0FVz2EQyIgJQCZj/3+6ePC3AslkLO20wrOM5qPoNoiI/DHwQeAC8AtB9UbgRMFuvUHdQvW9c9TP9533APcAdHR0FPcDomTVKnjzm00YRkaitsZZBiLzDFbNZGwa9ekWZH+rxZDi5pNyro1czjIJK3gcwmIs2k4Rke+JyL45yh4AVf2Eqm4CvgF8LGyDg+98QFU7VbWztbW1HF8ZHitWmDBMR0OdyqW+3s5loUuhtRV+/ucT/RBJDdPnN+HnctGegqq+Y4n/6xtYLOB+oA/YVLCtPajrw1xIhfU/COrb59g/HeRycMcd8OKLcPRo1NY410J7u/mY5woA1NbCnXfCoUPw8suVnfKUVlpbLXOwAucyWi5FeTRFpDDsvgd4KXj9KPBBMe4ALqjqKeBx4C4RaRaRZuAu4PFg27CI3BFkLn0QeKQY2yoOEbj5ZvNVelpJ5ZDL2cPida9bOCIsYqmLb3pTxc6Jk1p27KjYye2uhWJjCn8iItdjKak9WOYRWI/h3UA3lpL6YQBVHQzSWJ8N9vvkdNAZ+C1mUlK/S5KDzAvR3m6jj/bu9ThD3Fm71hZVWs5Mcs3N5k7avx+OHQvNNKcE5HIm9m1tUVtSVqTShwJ0dnZqV1dX1GaUnqkpeOEFOHFi8X2d8pLLWTB506bF912Ic+fg+edtoIQTL1paTBASPCxcRPaqaufses+AjivZrK3r2tZmDw7PXokHGzaYm6+mpvj/tWaN9RpefhkOH/ZYQxzIZGx08vbtqR0e7qIQd9avt6ktn3/eliFzoqG+3gLJpc52y2ZtIrWNG61nODi4+GeccFi1yhpijY1RWxIpLgqVQE2N5UafPAn79vlgt3JSVWXTGGzbFu5Is8ZGC0L39Vm8wVcRKh/ZrCUBpLh3UIiLQiWxYYO1VA8etCCluxvCQ8RiBjfcUBpX0VLZuNEmTuzuNpdSqebTcOamrc3cgbPnOkkxLgqVRi5nF3FHh41rOHcuaouSx7p1JgYrV0bz/dms+bU3b7Z4w/Hj3gAoNQ0Nliywdm3UlsQOF4VKpbHRBkSdOQMHDtjq9k5xtLbaw7i5OWpLjNpaeO1rza1x6BD09ro4FEtNjbmKOjp84ql5cFGodNrarLXT12etSk9vXD5tbRY3iIsYzKa+3gKgO3eaW6m316ddXy65nInr1q0+7ewi+NFJAiI26G3jRhOHQ4d84NtiZDJ2vLZvj85NtFzq622w3PXX23QoPT2+KPVi1NRYksCWLS4GS8SPUpKYFof2dktfPXIEBgaitipe1Naar37z5vIGkEtJba2lsV53nfUajh519+FsGhpMDNrbfdqYZeKikFTWrbMyPGyZSn196V3QR8RcbB0d5ipKStphNjsjcIODFpA+eTK9GUsidn63bCn9eJIU4aKQdBobLVi5axecOmXTZqSl97Bq1YxbrVJ7BUtl9WorN99s57mvzzLT0hCYXrnSzvOmTck/z2XARSEtVFXZTbNpE1y6ZC3KkydhaChqy0pLU5ONAl+/Pp2zkRae57ExcyOeOmUNgSQFp+vr7Rxv3Jj6EcilxkUhjdTVWYB1+3YbOXvmjJWBgcpzMeVy5ipobTXXgbcUZ6ipmXEvTUxYz+HsWejvt4ZBJSFiPaG1a+08V0pyQAXiopB2CgOv+bz1HAYGrJw/Hz+RqK211NHVq20my8bG5MQIwiSXm+lBgaUuF57nuKUyZzLm/lu92iYOXL3as4fKhB9lZ4ZMZsY3vTNYP2lkxIRieNjKyEh5WpmZjLkIVq60smqVleWsXeDMT329lek1zicm4MIFK8PDls108WJ5GgXV1ZYt1NhoZdUq++uDyyLBRcFZmIYGK4VMTcHoqJVLl8wFNTZm03tPTFiZmrKeRz5vwU4RK5mMZc1UVVmprrZSU2Nurbo6m4emttZ7AOUkl7MW+ey1wi9fvvo8T5/r6fM8OTlzrlVffa6nz3MuN3Oua2tnznV9fWpWNKsUSiIKIvJ7wP8AWlX1XLCk5uex1ddGgQ+p6nPBvncDfxR89NOq+mBQ/3pmVl57DPi4VvoKQEklm51pwTvJprbWe2cpo+j+mYhswtZaPl5Q/S5gZ1DuAb4c7LsauB94I3A7cH+wVjPBPr9e8LndxdrmOI7jLI9SOO0+B/w+UNiq3wN8XY2ngCYRWQ+8E3hCVQdV9TzwBLA72Naoqk8FvYOvA+8tgW2O4zjOMihKFERkD9Cnqj+btWkjULi4cG9Qt1B97xz1juM4ThlZNKYgIt8D1s2x6RPAH2Kuo7IiIvdgbik6prMnHMdxnKJZVBRU9R1z1YvIa4CtwM8srkw78JyI3A70AZsKdm8P6vqAt86q/0FQ3z7H/vPZ9ADwAEBnZ6cHox3HcUrENbuPVPUFVV2rqltUdQvm8rlNVU8DjwIfFOMO4IKqngIeB+4SkeYgwHwX8HiwbVhE7ggylz4IPFLkb3Mcx3GWSVjjFB7D0lG7sZTUDwOo6qCIfAp4Ntjvk6o6GLz+LWZSUr8bFMdxHKeMSKUPBejs7NSurq6ozXAcx6koRGSvqnbOrvdx5I7jOM4VKr6nICL9QE8EX70GOBfB9xZDJdoMbne5cbvLS1R2b1bVV61GVPGiEBUi0jVX1yvOVKLN4HaXG7e7vMTNbncfOY7jOFdwUXAcx3Gu4KJw7TwQtQHXQCXaDG53uXG7y0us7PaYguM4jnMF7yk4juM4V3BRcBzHca7gojAPIvLbIvKSiLwoIn9aUH+fiHSLyEEReWdB/e6grltE7i2o3yoiTwf13xSR6jLY/nsioiKyJngvIvKFwIbnReS2gn3vFpFDQbm7oP71IvJC8JkvBHNShWXvZ4Nj/byIfFtEmgq2xf54z8V89kVkyyYReVJE9gfX88eD+tUi8kRw7p+YXvDqWq6XkO3PishPROQ7wfs5z7GI1ATvu4PtWwr+x5zXUYg2N4nIt4Lr+oCI3FkpxxtV9TKrAL8AfA+oCd6vDf7uAn4G1GAzxB4GskE5DGwDqoN9dgWfeRh4f/D6K8Bvhmz7JmziwR5gTVD3bmwuKQHuAJ4O6lcDR4K/zcHr5mDbM8G+Enz2XSHafBdQFbz+DPCZSjne8/yeee2L6Hpej01WCbASeDk4tn8K3BvU31tw3Jd9vYRs/+8Cfwt8Z6FzjM2f9pXg9fuBby50HYVs84PAfwpeVwNNFXO8o7pQ41yCi+4dc9TfB9xX8P5x4M6gPD57v+Aknyt44F21X0i2fwu4BTjGjCj8JfCBgn0OBg+KDwB/WVD/l0HdeuClgvqr9gvZ/l8BvlEpx3ue3zCnfeW2YwH7HgF+cfo6COrWAwev5XoJ2dZ24PvA24DvLHSOp6+P4HVVsJ/Mdx2FaPMq4ChBIs/s4xjn462q7j6ah+uAtwRd0B+KyBuC+uWuKNcCDKnq5Kz6UJBkrIT3a8zMkBvr470A89kXOYFL5XXA00Cb2rT1AKeBtuD1co97mPw5ttxvPni/0Dm+Yl+w/UKwf7nt3gr0A/8zcHv9lYjUUxnHO7Sps2OPLLyiXBXWZbsDeAPwsIhsK6N587KI3ZGshLcUFrJbVR8J9vkEMAl8o5y2pQURaQD+HvgdVR0uDBOpqopIrPLTReQ9wFlV3Ssib43anmVQBdwG/LaqPi0in8fcRVeI4/GeJrWioPOsKAcgIr8J/INan+0ZEcljk1bNt6Ic89QPAE0iUhW0XBZcUa4YuyXClfCKsbvA/g8B7wHeHhx3iMHxvkYWsjsSRCSHCcI3VPUfguozIrJeVU+JyHrgbFC/3OslLN4E/LKIvBuoBRqBzzP/OZ62u1dEqjA3zgDlPx+9QK+qPh28/xYmCnE/3kbY/qlKLMBvYAsAgbmSTmC+yZu4OmB1BAsqVgWvtzITWLwp+PzfcXVQ7LfK9BuOMRNT+CWuDmQ9E9SvxnyfzUE5CqwOts0ONL87RFt3A/uB1ln1FXO8Z9k9r30RXc8CfB3481n1n+XqwOefXuv1Uobf8FZmAs1znmPgo1wdaH54oesoZHv/Bbg+eP3fgmNdEcc7kos07iW4kf8XsA94DnhbwbZPYNkLBynIyMEyCF4Otn2ioH5b8IDtDi7mmjL9hmPMiIIAXwxsewHoLNjv1wLbuoEPF9R3Br//MPAXzAqaldjWbkx4fxqUr1Ta8Z7jN81pX0S2vBlQ4PmCY/xuzN/+feAQlm033SBY9vVSht/wVmZEYc5zjPUm/i6ofwbYtth1FKK9twJdwTH/P8FDvSKOt09z4TiO41zBs48cx3GcK7goOI7jOFdwUXAcx3Gu4KLgOI7jXMFFwXEcx7mCi4LjOI5zBRcFx3Ec5wr/H+I7OId1+HeFAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "fig, ax = plt.subplots()\n", "# for maximally rotating black hole\n", "ax.fill(Xh1, Yh1, 'b', Xe1, Ye1, 'r', alpha=0.3)\n", "ax.fill(-1*Xh1, Yh1, 'b', -1*Xe1, Ye1, 'r', alpha=0.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot for rotating(normally) black hole" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[,\n", " ]" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD4CAYAAAAD6PrjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3de5BcZ3nn8e8z3XMfjeaikTS6S5acIGMMRhgIycZcFgyhYrYqtQW1tXESdl2bkFRSSRXB4Q+2kq2tJalaEjYXQoXswi5b4JBkoSiyrLkFcsFGNka25YtGsmyN7jMjzf3WPe/+8ZyeaY3n0udM3/v3qTo13afPmXlb3TrPed/nvVgIAREREYCmShdARESqh4KCiIgsU1AQEZFlCgoiIrJMQUFERJalK12ArdqxY0c4dOhQpYshIlJTHn/88ZEQwsDq/TUfFA4dOsTJkycrXQwRkZpiZi+ttV/NRyIiskxBQURElikoiIjIMgUFERFZpqAgIiLLFBRERGSZgoKIiCxTUBARkWU1P3hNpGhCgMVFWFggO7fI/NQii7MZMnMZFmczZOczZOazZBeyLC1GW2aJkF1iKRsIS74RrVGyhNGEP7Ym8y3VRFPKsHSKVHMTTc0pUi0rW7otTXNHM83taVo6m2nuaIbmZmhpgSbdw0npKShI/ctmYW4OZmdhbo7M1BxzN+eYG59nYWKOhakFMtPzLEwvsrjocWFpqQh/1tKkQmZLv8NsJSak29I0d7bQsq2Vlm2ttHa30tbTRkt3G7S1QXu7/2xu3nrhpWEpKEh9mJuD6WnfZmYIU9PMjc0wOzbL3Pg8c3N+yPy8X/RrRQiwsOAbUxkYyQAztxyTSkFrq8eE1lZo726mrbed9v4Omro6oKMDurqgs9MPMqvIe5HaoKAgtWVmBiYnV7apKbLjU8xMZJiegdkZP2R2tjh3+7Ugm/X3PJOLFZcWgUWwCdrbPCa0t3tM6Ohqormn04PEtm0rW2enmqcEUFCQarW05Bf98XHfJiZgYoKwmGF2FqamfJuehtk5QEuNv1Lw4Dg7m79zidbWSTo7J+nqukxXlweNpnSTB4rt26G7G3p6/Gdal4hGo09cqsP0NNy8CTdu+DYxAUtLLC35S3kVA7LZShe2ts3P+zY25s/NoLNzie7uCbZtm6Cry5ukAK9F9PRAb69v27ap+anOKShI+YXgF/3RUd9u3PCrVGRmZqVyMDnZOM1AlRLCSs0LAIOuTq8obN8+Sef4JE0XLvhr6bQHh74+33p78yKI1AMFBSmPiQkYGfFtdBQyK71ystmVVqKbN2srEVyX8oLEpUt+zc+1KG3fnqElcx2uX/djm5o8MPT3w8CAH6TcRE1TUJDSWFiAa9f84nH9+i01gdzLuZaiycnlrv1ShbLZlc8KPCfd2ws9vdDRvrRS43vhBa9J9PfDzp2+dXRUtvASm4KCFM/EBFy96tvNm6+40i8seDv22FheU4XUnFzP3+FhHxaRa0nq6MBrgLnvAHjyeudO2L3bD1I+ouopKEhyIfjt4+XLcOVKXp/IFZmMB4GREQWCejQ3501Mly55t9f+ft9aW6MDcu1Q5875oLpdu2Bw0AOFmpmqkoKCxDc2BhcvejBY1SwEHitu3oTrIzD+ygqD1KnZWa89DA97BWFgwCsHy3noxcWVA9JpDxB79ihAVBkFBSnMxIT/Z750aXXH92Vzc54+GBlRsrjR5SoIL73kgWFgwHuzLstk/Mbi4kWvQQwOwt69Xs1QE1NFKSjI+hYWPBBcuOBBYQ25FqRr19Y9RBrY0tJKp7OODg8OO3as6sW6uAgvv+xbezvs2wf793tGW8pOQUFuFYLf7r/8sicL1xkksLjoh129qlqBFGZmxmsOw8MeGHbt8kT1LWZn4cwZ3/r64OBBr0VoLETZKCiIm59fuVtbI2GcMzfnqYTRUQ0qk2Sy2ZUOSj29MLh7VdNSTq6r2tNPe83h0CHVHspAQaHRjY97z5BLlza8yk9NeTC4cRPNMyRFc/OGb11dXiHo6VkjpbC46N/Rc+e8/enIEf+p3ENJKCg0ohD8Nu3s2ZUJcNYxOem5QOULpJSmprzFqKPDg8O6QxpygyE7O+G22zz/oKalolJQaCRLS96gOzTko482oGAglTAz4/cqly55b9V1g8P0NJw6Bc89B4cP+6bFhYpCQaERZLOe4Tt71pMCG8iNVB0fL1PZRNYwO+tf18uXYe8+6O1Z58CFBXj+eT/44EGvPSyPnJMkFBTqWTYL58/7f5g1Bpnlm5/3YDA6Wp6iiRRiZgbOvOA5hwMH/OeaMhn/nr/4oiekjx5VcEhIQaEeLS15zeDMmU2DQSbjVfWrVzXyWKrX1BScPu3NSfv3b3C9X1ryhPRLL3mT0tGjalaKqWhjy80sZWY/NLOvRs8Pm9mjZjZkZl80s5Zof2v0fCh6/VDe73go2v+8mb2rWGVrGCH47f63v+3d+DYICLnhCKdO+bRFCghSC8bG/Ds7PLzJYkvZrOfOvvlNvznSykwFK+aEI78OPJv3/OPAJ0IIR4EbwAej/R8EbkT7PxEdh5kdB94P3AHcB/ypmalbQaFGRuC734Uf/nDDcQawctf14ou3LGsgUhNC8NrtU09t2nnOu7M+9xx861s+Ml93P5sqSlAws33AzwB/ET034G3Al6JDPgu8L3p8f/Sc6PW3R8ffD3whhDAfQngRGALuKUb56tr0NPzgB/DP/7xpV6FMxlMMp5/dtPORSNVbWPDKwPPPb9pK6h0snnwSvve9AiJJYytWTuEPgQ8DuXGJ/cDNEELuPnQY2Bs93gtcAAghZMxsPDp+L/D9vN+Zf84tzOxB4EGAAwcOFOkt1Jhs1qvFZ88WNLR4bMybWTUlhdSb8XGvNezd68s2bDimbXwc/vEf/eDjx9eYZ0O2XFMws/cC10IIjxehPAUJIXw6hHAihHBiYGCgXH+2ely54nmDM2c2DQiLi37Y0JACgtSvpSVvHTp9GmbWnsT3Vhcv+v+hc+fUpLRKMWoKbwF+1szeA7QB3cAfAT1mlo5qC/uAi9HxF4H9wLCZpYHtwGje/pz8cwS8CvzUUx4UCjA66rUD5Q2kUUxPwzNPe0VgcHCTWkMmA88841nru+6C7dvLVs5qtuWaQgjhoRDCvhDCITxR/K0Qwr8Bvg38XHTYA8CXo8dfiZ4Tvf6tEEKI9r8/6p10GDgGPLbV8tWNl17yO5sCAkKuy/bZswoI0nhynfCefbaAXAN4k9L3vucnqJdSSccp/DbwBTP7T8APgc9E+z8D/E8zGwLG8EBCCOEZM3sYOA1kgA+FEPQJzc56gmxkpKDDJyc9GCwslLhcIlVuasp7Zh886FN1bygEb2O9cgVe+1ro7S1LGauRhRpvTztx4kQ4efJkpYtRGhcu+Le6gNv9EHxKgOGLaBbTKpG1NKmgqlo16O/3gc4FzZ1n5oPebr+9rpcJNbPHQwgnVu/XiOZqtLgIP/qRX+ULPPzsWU1eJ7Ke0VHPNxw9Bh3tmxwcgvfOuH4d7r674dZwqN8wWKvGxuDv/77ggDA15bkyBQSRjc3NwelnYszvdfOm/1+8cKGk5ao2qilUk6EhH31ZYJPe9es+GK3GWwBFymZpyWvVU9NwYH8B6/Rks57TGx2FO+9siLUbFBSqweKiT09x9WpBh4ewsoSyiMR39QrMznjqIF3IVfDCBa85vOENdd+cpOajSpuc9O5wBV7hs1l44QUFBJGtmpjwwW6bLDGyYnLS5xer8/98CgqVdOUK/MM/FDwR0cKCd6XWAjgixTE354FhaqrAEzIZeOwxT0TXKQWFShka8onsChxdNjPrCeVNJkAVkZgyGU/l3bgR46TnnoMnniho3rFao6BQbiF4d9Nnn9382MjkJDx7WnMXiZTK0hKcGfLOGwW7eNFnJ66zkaIKCuWUycCjj3qWuEDj4z41sEbfi5RY8DVGCpxazI2N+ayrdVSFV1Aol4UFv6uIcSty46YnleuwhipStV5+2RfxKdjUlAeGOhkspKBQDnNz/qW5ebPgU27cgKEzGoMgUgnDw946VLC5Ofinf4qZmKhOCgqlNjPjAaHg7g1eQxgaUkAQqaSLF2PWGBYXvTWg4CHT1UlBoZRmZvzuIUZ74/i4aggi1WJ4uOAZZ1w263nDAmc1rkYKCqWSCwizhSwD5SYnvfuzAoJI9bhwIWavpGzWxzLUaI1BQaEU5ua8GhkjIMzMKqksUq1ePB8zXZCrMdRgjkFBodhyvYxiNBktLMAL6nYqUr2iNXgmJ2OckwsMNdYrSUGhmHLVxhhJ5dxcRnU2/kWk7uSWWSh4riTw5PP3v19T4xgUFIolBHj88VjVxRDg7Lma+r6INLRMxm/iYq19Pj/vNYYaufNTUCiWp5+OPXvixYtws/aaHEUa2txcgi7jU1M+11kNJA0VFIrh/HnfYhgbi9kHWkSqxsSEd1eNZWwMTp0qSXmKSUFhq0ZGvJYQw+ysz7EiIrXr8uUEnYsuXIBz50pSnmJRUNiK2VnPI8SoR2azXvVUTyOR2nfunKcMYjl9uqoHtykoJLW0BCdPxk4evfxyrOELIlLFcjd5sVIFIfhaDLG6MZWPgkJSzz4ba4I78KpmrJGRIlL1pqcT5Afn5z0wVOH0BQoKSVy9GrtdcHFReQSRenXpcsyBbeDTYFThsp4KCnHNz8OTT8Y+7cXzMfs2i0jtCH6fGDtX+MILVTcVhoJCXKdOxc4jjI5qPIJIvZufh+E4azCANx/98IdV1fNEQSGO4eGYa/V5s9FLL5WoPCJSVa5ejTXLjZuejrVme6ltOSiY2X4z+7aZnTazZ8zs16P9fWb2iJmdiX72RvvNzD5pZkNmdsrM7s77XQ9Ex58xswe2Wraimp+PPR4BPI6o2UikQQQfxxo7f/ziiz64rQoUo6aQAX4rhHAceBPwITM7DnwE+GYI4Rjwzeg5wLuBY9H2IPBn4EEE+BjwRuAe4GO5QFIVnnnGb/tjmJpSbyORRjMzA9euJTjx1KmqmAZjy0EhhHA5hPBE9HgSeBbYC9wPfDY67LPA+6LH9wOfC+77QI+ZDQLvAh4JIYyFEG4AjwD3bbV8RTEyEnPBVqdmI5HGNDwc+x7Suy9VwWjnouYUzOwQ8DrgUWBXCCG3kN0VYFf0eC9wIe+04WjfevvX+jsPmtlJMzt5vdS34iEkajYaGfGmQhFpPNmsd1ON7YUXKj6orWhBwcy6gL8GfiOEcMuqEiGEABRtlEYI4dMhhBMhhBMDAwPF+rVrO38+dgfkpaUEk2WJSF25ejXB9T2bheeeK0l5ClWUoGBmzXhA+HwI4W+i3VejZiGin7lWtovA/rzT90X71ttfObnJ02O6dq1mpk4XkVIJiVqdfdK8Cq7WVozeRwZ8Bng2hPBf8176CpDrQfQA8OW8/T8f9UJ6EzAeNTN9HXinmfVGCeZ3RvsqZ2go9tU9m9WU2CLiRkd9/fXYKthFNV2E3/EW4N8CT5lZbqjv7wD/BXjYzD4IvAT86+i1rwHvAYaAGeAXAUIIY2b2e8APouN+N4RQuT5aCwuJ5qW4dk1dUEVkxaWLcPRozJOuXfMuqn19JSnTRrYcFEII/wDYOi+/fY3jA/ChdX7XXwJ/udUyFcXZs7Gv7ktLsce2iUidGxvz3EJbW8wTn38e3vzmkpRpIxrRvJbFxdgrqYH3OIrdDU1E6t7lJD2RRkYqMi+SgsJazp9P1AakWoKIrCXxDePQUNHLshkFhdWWlhLlEsbHK969WESqVAgJZze4cqXsA54UFFa7dCnB+npwNcmwdhFpGNeuJVxTp8wLsSgorJYgl7CwEHsRNhFpMAsL3qIQ24ULZZ1aW0Eh38REosTOyAhFHK8tIvUqURNSJpNwFFwyCgr5LlzY/Jg1jIwUuRwiUpdu3kw4jinhtSkJBYWckGxM+syMEswiUpgQEvYyHRvzi00ZKCjkjIwkSjBXyboYIlIjRkcTnlimJiQFhZyEExYpKIhIHBOTCZuQyjSpmoICeJ0uwcizuTk1HYlITCFhL6SJibI0ISkogDfyJZjrWt1QRSSJxNeOMkyboKAACRdUTRjtRaThjY8nHMiW8FoVh4ICJOo8vLQUe0E2ERHAcwqzSdZZGB31i08JKShkMolu+WdmSv7ZiEgdm0hyU7m0VPKZUxUUxsYS1eNUSxCRrZhKeg0pcZdHBYWEUXdqqsjlEJGGknjyUwWFEkvYDaBMgwtFpE7NzydcY6HE3R4VFBLkEzKZRIOfRURukSjZvLCQ8MTCNHZQmJ9PdHUv4echIg0kcYvDxERRy5GvsYNCwmyxRjGLSDEkvsEsYVKzsYNCwkyPgoKIFEPiZmgFhRJJWHdTPkFEiiHxtaSEPV0UFBJIME2SiMgrJL6WKNFcIgnbgRQURKQYQkjYLbWEbdiNHRQS1t0SzYUuIrKGxSTXk2y2ZBeixg4KCW75l5Y055GIFE826bW9RE0WjR0UEkTabLYE5RCRhpVJek1plJqCmd1nZs+b2ZCZfaRkf2hpKdFEeAoKIlJMSwoK6zOzFPAnwLuB48AHzOx4Sf5Ywqv7UpKFMURE1pFosR0o2R1qVQUF4B5gKIRwLoSwAHwBuL8kfynpJ6GgICJFlDgoJD5xY9UWFPYCF/KeD0f7bmFmD5rZSTM7eT3BqmkiIrK2agsKBQkhfDqEcCKEcGJgYCDZLzFLeF6y00REiirpNWwT1RYULgL7857vi/YVX1Oyt96koCAiRZTwUgSpVFHLkVNtQeEHwDEzO2xmLcD7ga+U5C8l/AdN/AGKiKyhKem1vURBIV2S35pQCCFjZr8KfB1IAX8ZQnimZH+wuTn2GPMSfQ4i0qBSSW8006W5fFdVUAAIIXwN+FpZ/ljCoGBWssS/iDSYxNf2lpailiOnsRtDWlsTnVaiAC0iDai5OcFJTU0JTyzgV5fkt9aKtrZEp5UoQItIo7GEN5kJb2gL0dhBob090WnNCgoiUgQtzQl7lia8dhWisYNCR0ei00oYpEWkgSS+liS8dhWisYNCZ2ei09oUFESkCBK2YENXV1HLka+xg8K2bYlOK2HNTUQaiIJCtWlvT5TBV1AQkWJIfC3p7i5qOfI1dlCARP+4zc0l6w0mIg0kUWognVZOoaR6ehKdljAdISIC+I1lou7tPT0lmwwPFBSgtzfRaQoKIrIVnUnTAgmvWYVSUOjrS3Rawhy1iAgAXUlvLBNeswqloNDamugK39lZ0hqciNS5RLliMwWFskiwUE8qpSYkEUkmlUqYK+7tLfnkawoKADt3Jjpt+/Yil0NEGsK27oRrsyS8VsWhoADQ358o+ibsuCQiDa4n6Q3lrl1FLcdaFBTAQ3aCCNzZqfEKIhJfohvKjo6SDlrLUVDIGRxMdFqJe4eJSJ3p6ko4PmHPnqKXZS0KCjm7diVqQipxRwARqTOJrxl79xa1HOtRUMhJpRLVFrZtUxOSiBQuUVDYtq0sTUegoHCr/ftjn2IGO3aUoCwiUne6uxM2HR04UPSyrEdBIV9/f6IpaRUURKQQCYZEeUeYffuKXpZ1/1zZ/lKtOHQo9int7Zr2QkQ2lkol7JiyZ09ZF4ZXUFht//5ECecyjCkRkRo2MJBwwNrhw0Uvy0YUFFZLpxO13/X2KuEsIuuwhOPO+vrKPkpWQWEtR47EDulNTWUZbCgiNai3x+fejO3o0aKXZTMKCmtpb0+U2Nm509sNRUTyJRob291dkTtNBYX1HDsWe27sdFq5BRG5VXd3ok6N8GM/VvSyFEJBYT0dHYlyC7t3J0wmiUhdSjQQeft2v5hUwJYuX2b2B2b2nJmdMrO/NbOevNceMrMhM3vezN6Vt/++aN+QmX0kb/9hM3s02v9FMytfH6z13H577Pag5uaKfZYiUmW6uxN2V3/Vq4pelkJt9Z72EeDVIYTXAC8ADwGY2XHg/cAdwH3An5pZysxSwJ8A7waOAx+IjgX4OPCJEMJR4AbwwS2Wbeva2jzpHNPu3cotiEjCMWc7dyYc5VYcWwoKIYT/F0LIRE+/D+T+Ce4HvhBCmA8hvAgMAfdE21AI4VwIYQH4AnC/mRnwNuBL0fmfBd63lbIVzdGjHhxiSKfLNqGhiFSpvr4EuQQzOH588+NKqJit378E/F30eC9wIe+14Wjfevv7gZt5ASa3f01m9qCZnTSzk9evXy9S8deRTif6kHbtStgFTURqnlmiqdR8oFqFp0fYNCiY2TfM7Ok1tvvzjvkokAE+X8rC5oQQPh1COBFCODFQjmrW3r2xJzhqairrHFYiUkUGBxPcFLa1VazHUb5N53MIIbxjo9fN7BeA9wJvDyGEaPdFID9O7ov2sc7+UaDHzNJRbSH/+OrwmtfAd74DS0sFn9Lb650IxsdLVywRqS6trQnHJbz61Ymm2Cm2rfY+ug/4MPCzIYSZvJe+ArzfzFrN7DBwDHgM+AFwLOpp1IIno78SBZNvAz8Xnf8A8OWtlK3oOjsTRfFDh9RFVaSRHDiYoKPJ4GDi1R+LbauXqz8GtgGPmNmTZvYpgBDCM8DDwGng/wIfCiFko1rArwJfB54FHo6OBfht4DfNbAjPMXxmi2Urvttuiz0PSWtr2RZMEpEK6+vzKS1iaWmBO+8sSXmSsJUWn9p04sSJcPLkyfL9wakp+O53IZst+JQQ4PRpmJ4uYbmk6mQtTWq574TUu3Tar+2xJ8Y8caIitQQzezyEcGL1fjVsxNXVFbs3kpkPd4g5a4aI1JBDhxIEhAMHqqbZKEdBIYlDh2IPW25vh/3qjSRSl/r7E6y93NXlyeUqo6CQ1Gtf6/MjxbB7l/dGEpH60dqaYMHGVApe//qqnPpAQSGp5mZvC4zZtejIES3GI1IvzLz/Sexr+2te4xMjVSEFha3Yvh3uuivWKc3N/iVC+QWRmrf/QIKpLA4fTjgpUnkoKGzVvn2xJ83r7ob91fudEJEC9Pd7k3AsO3ZUfG6jzSgoFMPx47FXSBoc9C+ViNSejo4EeYTOzkRNzuVW3aWrFWZw992xs8iHDsXOVYtIhaXTvjBjrDxCSwu88Y01kVBUUCiWdNo/9BhX+VTK1/Gpge+JiOD3f8eOxZzsLpWCe+7xmkINUFAoptZWeNObYn1jWlp8SqUq7JkmIqscORJzZmszbzLq7S1ZmYpNQaHYOjs9MMS4/e/o8LV8NOJZpHodOBAzD5hrVt65s2RlKgUFhVLo7vampBjT4G7fHnViUmAQqTqDgwnWXr/rrppcglFBoVR6e2MHhv5+OHSwhGUSkdgGBhKsonbXXQmXXqs8BYVS6uuL3ZS0c6dWbBOpFv39MbuemvkUODX8n1hBodR6e+HNb/aMcoF2767ZmwyRutHfH3N241wOocb/8yoolMP27fCWt/garAUaHKzpmw2RmhY7IOS6ndZgDmE1BYVy6eqCn/zJWP3Zdu9OMGpSRLZkYCBmQGhu9mbiGutltB4FhXJqb/caQ4x+bTt3agI9kXIZHPT56goOCB0dfrMXezGF6qWgUG65u4oY7Y79/XD7saqfMkWkpu3fHzMd0NvrASH2NKnVTZeZSmhq8h4Kr3pVwbckPT3w4z+uKTFEii23JkKsVTH37YOf+ImY813UBgWFSjp6FN7whoKv9LnlodvbS1wukQaRTvvNVsEtumb+n/B1r6vbqnt9vqtasmsX/NRPFZyAbm31CoaW9RTZmrY2v74X3PejpcWbfm+7raTlqjQFhWrQ2emBocDVmNJpn1019rB7EQG8Ofb48Ri9xPv64Kd/2hfJqXOFz8EgpZVKeZV0xw546inIZjc83MzHMXR2wosvwtJSmcopUuP27IG9e2P0MDp61KcyrtPmotUUFKrN/v3eq+GJJ2B8fNPD+/uhvQOGzsDcXBnKJ1KjUilv+enpKfCEtjbvEDIwUNJyVZvGCH21JjfQrcD5tDva4Y47tLynyHq6uuDVr44REAYHvbmowQICqKZQvZqaPKO8axc8+SRMT294eO4uqLsbXn5509YnkYaxZ49vBbX+NDd79Cgwv1ePFBSqXS7B9dxznjwIYcPDBwa8N8W5czA1VaYyilShlhafrqK7u8ATdu+GO++MNUdZPSpK85GZ/ZaZBTPbET03M/ukmQ2Z2Skzuzvv2AfM7Ey0PZC3//Vm9lR0zifNtA7ZslTK24fe8paC+s+1tXklI1YyTaSO9Pf7DX9BAaG1FV7/eh8z1OABAYoQFMxsP/BO4OW83e8GjkXbg8CfRcf2AR8D3gjcA3zMzHKLl/4Z8O/zzrtvq2WrO7298C/+hV/xN1nU2cyDwh13+PQsIo2guRmOHfOm1ILWtzp4EN761rqY3bRYilFT+ATwYSC/XeN+4HPBfR/oMbNB4F3AIyGEsRDCDeAR4L7ote4QwvdDCAH4HPC+IpSt/jQ1eQL6rW8taKBCR4f3x963T7UGqW8DA97609u7+bFs3+6dOV7zGs0ds8qWcgpmdj9wMYTwo1WtPXuBC3nPh6N9G+0fXmP/en/3QbwGwoFGXXSgvd2ru9evw9NPb5hAaGryG6G+Pjh/HiYmyldMkVJrbfWZTQtqKmpp8Zr2/v26S1rHpkHBzL4BrHVL+lHgd/Cmo7IKIXwa+DTAiRMnNs681ruBAbj3XnjpJXj+eVhYWPfQtjaf52VkBC5cgMXF8hVTpNhyNzu7dxfQs6ipybPOR4+qZrCJTYNCCOEda+03szuBw0CulrAPeMLM7gEuAvmT0O6L9l0E7l21/zvR/n1rHC+FMPPVePbtg6Eh73q0QZ/UHTu8ij18Ea5d3bRDk0jV6e31Ef0FTVK6b5/fDWkmyYIkbj4KITwFLC81ZGbngRMhhBEz+wrwq2b2BTypPB5CuGxmXwf+c15y+Z3AQyGEMTObMLM3AY8CPw/8t6Rla1i5KR8PHYIzZ3zAwjrzX6RScPAA7BzwWsPNm+UtqkgSHR0eDApqKtq1y/8/FNwnVaB04xS+BrwHGAJmgF8EiC7+vwf8IDrud0MIY9HjXwH+B9AO/F20SRJtbZ5xu+02Dw4XLqxbHWhv98n1JiY8hszMlLmsIgVobvYb/h07CkgF7NzpcxUVPHxZ8lmo8baDEydOhJMnT1a6GNVtZrN/uAsAAArBSURBVMablS5c2HTmvNFRGB6G+fkyla2OZS1NKmQqXYyalkp53mDnzk17YXvN4NixArsfiZk9HkI4sXq/RjQ3go4O73p3++1w9qxXCTJrX6z6+72X0sgIXLy4Yd5apGRSKU8g7969STAw86hx9KiaiYpEQaGRtLX5aLbbb/e+qefPrzm1qpl3aurv95rDpUuqOUh55ILBrl2bDD5LpTy5cOSIRmcWmYJCI8of9nnpks+ptEamuanJg8OOHR4cLl+G2dkKlFfqXnOzB4JduzapGXR0eEeKAwfUtbREFBQaWVOTZ+/27YMbN7zmcOnSK/IOZh4YduyAGzfhymWYnKxMkaW+tLX5LNX9/ZuMNRgY8BFqO3dq0FmJKSiI6+317dWv9oT0yy+veeXv7fFtehquXIGxMY1zkPi6u72ZaPv2Da7xra1eIzhwQE1EZaSgILdqbvZ22iNHvPbw8svebrRq+HNnp7c+HTjgM21cu6aktGwslfLa5s6dG4wja2ryAw4cUK2gQhQUZH35tYerV72v6vXrtzQvNTd754/BQbg5Dtev+U9Ue5BIV5e3/vT1bZAv6O31Zsw9e3x+IqkYBQXZXK6z+J49XmO4dMm30dHltiOzlaalhQXv0joyonWjG1Vzs+cJdgz4crFr6u72+d337FHzUBVRUJB4mpt9DvqDB72f6tWrKwEiqkG0tKzEkKkpf2lsTBPw1bumJr/h7+/fIFfQ0+PVysFBb4OUqqOgIMnlJwIzGU8sXLniP6MI0NXl24EDnrceHfVUxTpj56TGNDV5AOjr8+v9K5qHmpo8kZAbfKCVzaqegoIURzq9Uj0Iwa/81675Nj6OmbcWdHd7N/PJST/kxg0lqGtNKuUBoLfXA8IrAkFnpycRdu70gLDp/BRSTRQUpPjM/Naxr89nqcwlGa5fh9FRbHp6OUAcPOhTM9286dvUNEpSV6G2Ng8EPT1e87tlTEFrq7cZ5UY6Kj9Q0xQUpPTykwzgw6JHR5eTDR1M0dHhL2cyMD7us7ZOTGh6jUpJpz1ob9/uP29Zt6CtbWWSrP5+2LatYuWU4lNQkPJrb18ZSQ1ek4jaktI3btDfPk5/v+ck5udXAsTUlIJEqaTT0LUNurd5EGhvjxLFuaRBb69XE/r6tFhNnVNQkMpraVmZ+CZnehpu3qR1fJyBaGNxkYUFz0dMTcP0lB+mEdUxmXcT7exc6QjQ3o63/efa9Xp6PBhs21bAWpdSTxQUpDp1dvq2d+/Kvrk5WiYm6J+cpD+qOixNTDE7mWFmxgPEzIy3Tm2wGmlDMfMLfkcHdHRCZwd0dDWR6u70C35u6+72gzSCuOEpKEjtaGvzbefyKrA0AZ1zc3ROTTEwPe2RYXqaubEZZkdnmJ3MMDfH8lavXWGbmlb+edraoL2zifa+dtr6O2nq6lgJsl1duvjLhhQUpPblroQ7dqzsirbexcWV6sPcHIuTc8zdmGVufJ6FiTkWJudZnFlkYT6wsFC9QSOV8nGDLS3Q0tlMc2cLrdvbaN3eRtv2Vlp72qNo0O5ba6su/JKIgoLUt+Zmbxvfvt2fRtst/WVC8GT3wgJL84ssTC2wML3IwvQimdlFFmczZOczZBeyKz8XlwiZLEuZJUI22pYCYckTHGEpYKkmLNuENRmYYU2GpZpoSkU/m1Okmv1nU2uadEuKdFuadFua5o5m0h3NtHQ009rVTKq9JYoILWrjl5JSUBAx8zvr1laatkHbDq9liDQi3XKIiMgyBQUREVmmoCAiIssUFEREZJmCgoiILFNQEBGRZQoKIiKyTEFBRESWKSiIiMgyCzU+77CZXQdeKuKv3AGMFPH3VYLeQ3XQe6gOeg9rOxhCGFi9s+aDQrGZ2ckQwolKl2Mr9B6qg95DddB7iEfNRyIiskxBQURElikovNKnK12AItB7qA56D9VB7yEG5RRERGSZagoiIrJMQUFERJY1XFAws18zs+fM7Bkz+/28/Q+Z2ZCZPW9m78rbf1+0b8jMPpK3/7CZPRrt/6KZtZT5ffyWmQUz2xE9NzP7ZFSeU2Z2d96xD5jZmWh7IG//683sqeicT5qVZ1FfM/uD6DM4ZWZ/a2Y9ea/V1OewnvXKW2lmtt/Mvm1mp6P/A78e7e8zs0ei78gjZtYb7Y/9vSrje0mZ2Q/N7KvR8zW/C2bWGj0fil4/lPc71vy+lfE99JjZl6L/D8+a2Zsr/lmEEBpmA94KfANojZ7vjH4eB34EtAKHgbNAKtrOAkeAluiY49E5DwPvjx5/CvjlMr6P/cDX8UF7O6J97wH+DjDgTcCj0f4+4Fz0szd63Bu99lh0rEXnvrtM5X8nkI4efxz4eC1+Dhu8v3XLW+kNGATujh5vA16I/t1/H/hItP8jeZ9J7O9VGd/LbwL/G/jqRt8F4FeAT0WP3w98caPvW5nfw2eBfxc9bgF6Kv1ZVPxLWuYP4GHgHWvsfwh4KO/514E3R9vXVx8XfSgjeRe2W44rw/v4EnAXcJ6VoPDnwAfyjnk+ugB8APjzvP1/Hu0bBJ7L23/LcWV8L/8K+Hwtfg4bvKc1y1vpcq1T1i8D/zL3fYn2DQLPJ/lelbHc+4BvAm8DvrrRdyH3PYoep6PjbL3vWxnfw3bgRaIOP6v/jSv1WTRa89HtwE9FVci/N7M3RPv3AhfyjhuO9q23vx+4GULIrNpfcmZ2P3AxhPCjVS/FfQ97o8er95fbL+F3P1BDn8Mm1itvVYmaUV4HPArsCiFcjl66AuyKHsf9TMrlD4EPA0vR842+C8tljV4fj46v9Hs4DFwH/nvUDPYXZtZJhT+LdNITq5WZfQPYvcZLH8Xfbx9e9XoD8LCZHSlj8QqyyXv4Hbz5papt9B5CCF+OjvkokAE+X86yCZhZF/DXwG+EECby00khhGBmVdtX3czeC1wLITxuZvdWujxbkAbuBn4thPComf0R3ly0rBKfRd0FhRDCO9Z7zcx+Gfib4HWsx8xsCZ9o6iLeTp+zL9rHOvtHgR4zS0d3HvnHl+w9mNmd+N3Fj6L/xPuAJ8zsng3ew0Xg3lX7vxPt37fG8UWx0ecAYGa/ALwXeHv0eUCVfQ5bsNH7qDgza8YDwudDCH8T7b5qZoMhhMtmNghci/bH/V6Vw1uAnzWz9wBtQDfwR6z/Xci9h2EzS+PNNqNU/nMaBoZDCI9Gz7+EB4XKfhblaj+rhg34D8DvRo9vx6tcBtzBrQmnc3iyMB09PsxKwvCO6Py/4tak1q9U4P2cZyWn8DPcmoR6LNrfh7db9kbbi0Bf9NrqRPN7ylTu+4DTwMCq/TX5Oazx/tYtb6W36LP+HPCHq/b/AbcmN38/6feqzO/nXlYSzWt+F4APcWui+eGNvm9lLv/3gB+LHv/H6HOo6GdR8S9pmT+AFuB/AU8DTwBvy3vto3jvg+fJ64WDZ/xfiF77aN7+I9FFdSj6MrZW4P2cZyUoGPAnUTmfAk7kHfdLUTmHgF/M238i+rc4C/wxqxJeJSz3EB6Qn4y2T9Xy57DOe1yzvJXegJ8EAnAq79//PXgb+zeBM3gPvdyNQ+zvVZnfz72sBIU1vwt4beKvov2PAUc2+76VsfyvBU5Gn8f/iS7qFf0sNM2FiIgsa7TeRyIisgEFBRERWaagICIiyxQURERkmYKCiIgsU1AQEZFlCgoiIrLs/wMEaVICaBk/5AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "fig, ax = plt.subplots()\n", "ax.fill(Xh2, Yh2, 'b', Xe2, Ye2, 'r', alpha=0.3)\n", "ax.fill(-1*Xh2, Yh2, 'b', -1*Xe2, Ye2, 'r', alpha=0.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - The inner body represents event horizon and outer one represents ergosphere. It can be concluded that with decrease in angular momentum, radius of event horizon increases, and that of ergosphere decreases." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.6" } }, "nbformat": 4, "nbformat_minor": 2 }